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ABSTRACT 

Without imposing any theoretical models and assumptions, we present a multi-variable regression 
analysis to several observable quantities for a sample of 15 gamma-ray bursts (GRBs). The observ- 
ables used in the analysis includes the isotropic gamma-ray energy {E^^{so)i the peak energy of the 
vFi^ spectrum in the rest frame {E^), and the rest frame break time of the optical afterglow light 
curves(i|3). A strong dependence of E'^^iso on E^ and is derived, which reads i?-y^iso/10^^ergs = 
(0.85 ± 0.21) X (Ep/lOO keV)i-9'*±° " x (4/lday)-i-24±o.23 ^ Universe with VLm = 0.28 and 
Hq = 71.3 km s~^ Mpc~^. We also extend the analysis to the isotropic afterglow energies in the X-ray 
and the optical bands, respectively, and find that they are essentially not correlated with E^ and t^. 

Regarding the Ey^iso{Ep,ty^) relationship as a luminosity indicator, we explore the possible constraints 
on the cosmological parameters using the GRB sample. Since there is no low-redshift GRB sample 
to calibrate this relationship, we weigh the probability of using the relationship in each cosmology to 
serve as a standard candle by statistics, and then use this cosmology-weighed standard candle to 
evaluate cosmological parameters. Our results indicate 0.05 < Hm < 0.50 at la level, with the most 
possible value of flu being 0.28. The best value of is 0.64, but it is less constrained. Only a 
loose limit of f^A < 1.2 is obtained at la level. In the case of a flat Universe, the la constraints are 
0.13 < < 0.49 and 0.50 < < 0.85, respectively. The decelerating factor {q) and its cosmological 
evolution (dq/dz) are also investigated with an evolutionary form of g = go + zdq/dz. The best-fit values 
are {qo, dq/dz) ^ (-1.00,1.12), with -2.23 < go < 0.26 and -0.07 < dq/dz < 3.48 at la level. The in- 
ferred transition redshift between the deceleration and the acceleration phases is 0.78l;o'23 (Ic). Through 
Monte Carlo simulations, we find that the GRB sample satisfying our relationship observationally tends 
to be a soft and bright one, and that the contraints on the cosmological parameters can be much im- 
proved either by the enlargement of the sample size or by the increase of the observational precision. 
Although the sample may not expand significantly in the Swift era, a significant increase of the sample is 
expected in the long term future. Our similations indicate that with a sample of 50 GRBs satisfying our 
multi-variable standard candle, one can achieve a constraint to the cosmological parameters comparable 
to that derived from 157 Type la supernovae. Furthermore, the detections of a few high redshift GRBs 
satisfying the correlation could greatly tighten the constraints. Identifying high-z GRBs and measuring 
their E'p and ^j, are therefore essential for the GRB cosmology in the Swift era. 

Subject headings: cosmological parameters — cosmology: observations — gamma-rays: bursts 



1. INTRODUCTION 

Long gamma-ray bursts (GRBs) are originated from cos- 
mological distances (Metzger et al. 1997). Their births fol- 
low the star formation history of the universe (e.g. Totani 
1997; Paczynski 1998; Bromm k Loeb 2002; Lin et al. 
2004). GRBs therefore promise to serve as a new probe 
of cosmology and galaxy evolution (e.g., Djorgovski et al. 
2003). It is well known that Type la supernovae (SNe la) 
are a perfect standard candle to measure the local universe 
up to a redshift of ~ 2 (e.g., Riees et al. 2004). Gamma- 
ray photons (with energy from tens of keV to several MeV) 
from GRBs are almost immune to dust extinction. They 
should be detectable out to a very high redshift (Lamb 
& Reichart 2000; Ciardi & Loeb 2000; Gou et al. 2004). 
Hence, GRBs are potentially a more promising ruler than 
SNe la at higher redshifts. 

This issue has attracted much attention in GRB com- 



munity. Frail et al. (2001) found that the geometrically- 
corrected gamma-ray energy i?jet for long GRBs is nar- 
rowly clustered around 5 x 10^° ergs, suggesting that GRBs 
can be potentially a standard candle. A refined analysis 
by Bloom, Frail, & Kulkarni (2003a) suggests that i?jot is 
clustered at 1.3 x 10^^ ergs, but the dispersion of Ej^t is 
too large for the purpose of constraining cosmological pa- 
rameters. Schaefer (2003) considered two other luminos- 
ity indicators proposed earlier, i.e. the variability (Feni- 
more & Ramirez- Ruiz 2000; Reichart et al. 2001) and the 
spectral lag (Norris, Marani, & Bonnell 2000), for nine 
GRBs with known redshifts, and pose an upper limit of 
r^M < 0.35 {la) for a flat universe. Using 12 BeppoSAX 
bursts, Amati et al. (2002) found a relationship between 
the isotropic-equivalent energy radiated during the prompt 
phase (Ey^iso) and rest frame peak energy in the GRB spec- 
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trum (E'p), i.e., E'p oc £;^(4. This relation was confirmed 



1 



2 



and extended to X-ray flashes by HETE-2 observations 
(Sakamoto et al. 2004a; Lamb et al. 2005a). In addition, 
it also exists in the BATSE bursts (Lloyd et al. 2000), and 
even in different pulses within a single GRB (Liang, Dai & 
Wu 2004). Possible theoretical explanations of this corre- 
lation have been proposed (Zhang & Meszaros 2002a; Dai 
& Lu 2002; Yamazaki, loka, & Nakamura 2004; Eichler & 
Levinson 2004; Rees & Meszaros 2005). Because of a large 
dispersion, this relationship is not tight enough to serve a 
standard candle for precision cosmology, cither. 

Ghirlanda et al. (2004a) found a tighter correlation be- 
tween GRB jet energy and E^, which reads Ejet oc (Ep)^^'^, 
where Ejci = E^^isoi^ ^ cosOjci), and Ojct is the jet opening 
angle inferred from the "jet" break time imprinted in the 
light curves (usually in the optical band, and in some cases 
in the X-ray and the radio bands) by assuming a uniform 
top-hat jet configuration. It is puzzling from the theoreti- 
cal point of view how a global geometric quantity (jet an- 
gle) would conspire with £'^,,180 to aS'ect E'^. Nonetheless, 
the correlation has a very small scatter which is arguably 
fine enough to study cosmology. By assuming that the cor- 
relation is intrinsic, Dai. Liang, & Xu (2004) constrained 
the mass content of the universe to be Qm = 0.35 ±0.15 in 
the case of a flat Universe with a sample of 14 GRBs. They 
also constrained the dark matter equation-of-state param- 
eter in the range of u> = —0.84^0 ^3 at la level. Ghirlanda 
et al. (2004b) evaluated the goodness of this relationship 
in diflferent cosmologies by exploring the full cosmological 
parameter space and came up with similar conclusions. 
Priedmann & Bloom (2005) suggested that this relation- 
ship is only marginal but not adequate enough for a pre- 
cision cosmology study. The main criticisms arc related 
to several assumptions involved in the current Ghirlanda- 
relation, such as constant medium density (which could 
vary in different bursts, e.g. Panaitescu & Kumar 2002), 
constant radiative efficiency (which also varies from burst 
to burst, e.g. Lloyd-Ronning & Zhang 2004; Bloom et 
al. 2003a and references therein), and the assumption of 
the top-hat jet configTiration (in principle jets are possibly 
structured, Rossi ct al. 2002; Zhang & Meszaros 2002bj. 
Nonetheless, the Ghirlanda-relation has motivated much 
work on measuring cosmology with GRBs (e.g. Firmani 
et al. 2005; Qin et al. 2005; Xu et al. 2005; Xu 2005; 
Mortsell & SoUerman 2005). 

In this work, we further address the GRB standard can- 
dle problem by a new statistical approach. Instead of 
sticking to the jet model and searching for the correla- 
tion between E'^ and Ejot (which requires a model- and 
parameter-dependent jet angle), we start with pure ob- 
servable quantities to search for possible multi-variable 
correlations by using a regression method. Similar tech- 
nique was employed by Schacfcr (2003). The motivations 
of our analysis are two-fold. First, within the jet model, 
there is no confident interpretation to the Ghirlanda rela- 
tion. It is relatively easy to imagine possible correlations 
between E^ and -E^^iso (e.g. Zhang & Meszaros 2002a; 
Rees & Meszaros 2005), since the latter is also a mani- 
festation of the energy per solid angle along the line of 
sight, which could be possibly related to the emission spec- 
trum. However, it is hard to imagine how the global ge- 

^The simple Ghirlanda-relation could be derived from the standard 
is constant for all GRBs, which is not true (see also Wu, Dai, & Liang 



ometry of the emitter would infiuence the local emission 
property^. Since there is no straightforward explanation 
for the Ghirlanda relation, one does not have to stick to 
this theoretical framework, but should rather try to look 
for some empirical correlations instead. This would al- 
low more freedom for possible interpretations. Second, 
within various theoretical models (e.g. Table 1 of Zhang & 
Meszaros 2002a), the value of E'^ depends on multiple pa- 
rameters. The problem is intrinsically multi-dimensional. 
It is pertinent to search for multi-variable correlations 
rather than searching for correlations between two pa- 
rameters only. The Ghirlanda-relation is a relation that 
bridges the prompt emission and the afterglow phases. It 
is also worth checking whether or not there are similar 
relationships for other parameters. Below we perform a 
blind search for the possible multi-variable correlations 
among several essential observable quantities, including 
the isotropic gamma-ray energy E^^iso, the isotropic X-ray 
afterglow energy i?xA,iso, the isotropic optical afterglow 
energy E'oA.iso, the cosmological rest frame peak energy 
Ep, and the cosmological rest frame temporal break in the 

optical afterglow light curve {t^). We describe our sam- 
ple selection criteria and the data reduction method in 
§2. Results of multiple regression analysis are presented in 
§3. A strong dependence of E^^iso on E^ and fj, is derived 
from our multi- variable regression analysis. Regarding the 
E^^iao{Ep,t^^) relationship as a luminosity indicator, in §4 
we explore the possible constraints on cosmological param- 
eters using the GRB sample. In addition (§5), we perform 
Monte-Carlo simulations to investigate the characteristics 
of the GRB sample satisfying the relationship observation- 
ally, and examine how both the sample size and the ob- 
servational precision affect the constraints on cosmological 
parameters. Conclusions and discussion are presented in 
§6. Throughout the work the Hubble constant is adopted 
as Ho = 71.3 km s'^ Mpc'^ 



2. SAMPLE SELECTION AND DATA REDUCTION 

Our sample includes 15 bursts with measurements of 
the redshift z, the spectral peak energy Ep, and the op- 
tical break time tb. It has been suggested that the ob- 
served Amati-relationship and the Ghirlanda-relationship 
are likely due to some selection effects (Nakar & Piran 
2005; Band & Preece 2005). The sample from which the re- 
lations are drawn may therefore be ill-defined if the parent 
sample is the whole GRB population. However, we believe 
that due to the great diversity of GRBs and their after- 
glow observations, one does not have to require all GRBs 
to form a global sample to serve as a standard candle. If 
one can identify a subclass of GRBs to act as a standard 
candle (such as SNe la in the supernova zoo), such a sam- 
ple could give meaningful implications to cosmology. Our 
selected GRBs belong to such a category, which assemble 
a unique and homogeneous subclass. Since not all GRBs 
necessarily have an Ep or a, th, the parent sample of our 
small sample is also only a sub-class of the whole GRB 
population. Notice that in order to preserve homogene- 
ity, we do not include those bursts whose afterglow break 
times were observed in the radio band (GRB 970508, GRB 

afterglow model and the Amati-relation, but one has to assume that t'^ 
2004). 
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000418, GRB 020124) or in the X-ray band (GRB970828), 
but were not scon in the optical band. Since wc arc not 
sticking to the jet model, we do not automatically accept 
that there should be a temporal break as well in the op- 
tical band. We also exclude those bursts whose Ep or ti, 
are not directly measured (but with upper or lower limits 
inferred from theoretical modeling). This gives a sample 
of 15 bursts up to Feb, 2005. They are tabulated in Ta- 
ble 1 with the following headings: (1) the GRB name; 

(2) the redshift; the spectral fitting parameters including 

(3) the spectral peak energy Ep (with error aSp), (4) the 
low-energy spectral index a, and (5) the high-energy spec- 
tral index /3; (6) the 7-ray fluence (S^y) normalized to a 
standard band pass (1 — 10* keV in the cosmological rest 
frame) according to spectral fitting parameters (with er- 
ror {(TS-y))', (7) the corresponding observation energy band; 
and (8) the references for these observational data. Our 
GRB sample essentially resemble those used in Ghirlanda 
et al. (2004a), Dai et al. (2004), and Xu et al. (2005). 
These bursts are included in the Table 1 of Fricdmann & 
Bloom (2005), but that table also includes those bursts 
with only limits for Ep, tj, and z, as well as those bursts 
whose th were observed in the non-optical bands (or in- 
ferred from theoretical model fittings). We believe that 
our sample is more homogeneous than the sample listed in 
Friedmann & Bloom (2005). 

The X-ray and optical afterglow data of these GRBs are 
listed in Table 2 with the following headings: (1) the GRB 
name; (2) the X-ray afterglow temporal decay index, a^; 
(3) the epoch of x-ray afterglow observation (in units of 
hours); (4) the 2-10 keV X-ray flux {F^, in units of IQ-^^ 
erg cm s ^^) at the corresponding epoch; (5) the 2-10 
keV X-ray afterglow flux normalized to 10 hours after the 
burst trigger (including the error); (6) the temporal break 
(including the error) of the optical afterglow light curves 
(ib); (7) the optical temporal decay index before the break 
(ai); (8) the optical temporal decay index after the break 
(012); (9) the rofcroncos; (10) the R-band optical afterglow 
magnitudes at 11 hours. We find that the mean values of 
ax, oi\, and a2 in our sample are 1.41, 1.0, and 2.0, re- 
spectively. For those bursts whose a^, ai, and 02 values 
are not available, we take these means in our calculation. 

With the data collected in Tables 1 and 2, we calcu- 
late the total isotropic emission energies in the gamma-ray 
prompt phase (E^^iso), in the X-ray afterglow (i?xA,iso)) 
and in the optical afterglow (R band) (i^oA.iso), i-e. 



K ■ = 



A-KDl{z)Syk 
1+^ ' 



and 



l + z 



(1) 



(2) 



l + z 



. (3) 

Here Di^{z) is the luminosity distance at the redshift z, k 
is a k-correction factor to correct the observed gamma-ray 
fluence at an observed bandpass to a given bandpass in 
the cosmological rest frame (1 — 10* keV in this analysis), 
ti and t2 are, respectively, the starting and the ending 



times of the afterglow phase, is the flux of the X-ray 
afterglow in the 2-10 kcV band, and Fji is the flux of the 
optical afterglow in the R band {Fji = vr/v). Since the 
very early afterglows might be significantly different from 
the later afterglows, which were not directly detected for 
the GRBs in our sample, we thus take ti = 1 hour. We 
also choose t2 = 30 days. The derived E'^.iso) -E'xA.iso) and 
-EoA.iso arc tabulated in Tablc^ 3. 

3. MULTIPLE VARIABLE REGRESSION ANALYSIS 

As mentioned in §1, previous authors interpret the re- 
lationship among E'^^iso, -Bp, and based on the GRB 
jet model (e.g., Rhoads 1999; Sari et al. 1999). In this 
scenario, the relationship among the three quantities be- 
comes the Ejet oc (E'p)^/^ relationship. When this relation 
is expanded, one gets E~^_isot[^ cx (-Ep)^. The indices for 
Ej^iso and are not independent and are bound by the 
jet model. However, since the current jet model is difficult 
to accommodate the Ej^t — Ep relationship, we no longer 
need to assume an underlying correlation between -E-y.iso 
and t'-^. We therefore leave all the indices as free parame- 
ters and perform a multiple variable regression analysis to 
search a possible empirical relationship among E^,iso, Ep, 
and f^. We also extend our analysis to search the depen- 
dence of Ex.iso and Eraso on Ep and t^, respectively. The 
regression model we use reads 



lO'^°Ep''H^''^ 



(4) 



where E^ = Ep{l + z) and = ib/(l + ^)- We measure the 
significance level of the dependence of each variable on the 
model by the probability of a t-test (pt). The significance 
of the global regression is measuK^d by a F-test (the F- 
test statistics and the corresponding significance level pp) 
and a Spearman linear correlation between logEiso and 
logEiso (the correlation coefficient r and the significance 
level ps)- We find that Ej^iso strongly depends on both 
E and with a very small uncertainty (Table 4, Fig.l). 
The actual dependence format depends on the cosmology 
adopted. For a fiat Universe with Om = 0.28, this relation 
reads 



E^,iso,52 = (0.85 ± 0.21) X 



E^ 



1.94±0.17 



1 day 



100 keV 

1.24±0.23 



(5) 
(6) 



where -E-y,iso.52 = E-j,.iso/10'''^ergs. However, when we test 

the possible correlations among Ex,iso (or Eo.iso), Ep and 

t^, no significant correlation is found (Table 4, Figs. 2 and 
3). 

4. LUMINOSITY INDICATOR AND COSMOLOGICAL 
IMPLICATIONS 

The dispersion of the E-ySso{Ep,tfj) relationship is so 
small that it could potentially serve as a luminosity indica- 
tor for the cosmological study. This relationship is purely 
empirical, exclusively using directly measured quantities, 
and without imposing any theoretical models and assump- 
tions. It therefore suffers less uncertainties/criticisms than 
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docs the Ghirlanda relation (e.g. Friedmann & Bloom 
2005). Below we will discuss the eosmological implications 
for this new empirical luminosity indicator. 

The distance modulus of a GRB, which is defined as 
/i = 51og(DL/10pc), could be measured by this luminos- 
ity indicator as 

fi = 2.5[Ko + Ki\ogEp + K2^ogti, (7) 
- log(47rS'^fc) + log(l + z)] - 97.45. (8) 

Since the luminosity indicator is cosmology-dependent, /t 
is also cosmology-dependent. We therefore cannot directly 
use this relationship for our purpose. Ideally, it should be 
calibrated by local GRBs (e.g., z < 0.1), as is the case of 
Type la supernova cosmology. However, the GRB low red- 
shift sample is small. More importantly, the local GRBs 
appear to have different characteristics than the eosmolog- 
ical ones (e.g. long lag, less luminous etc), so that they 
may not belong to the subclass of GRBs we are discussing. 
We are left out without a real (cosmological-independent) 
luminosity indicator at this time. 

We adopt the following method to circumvent the dif- 
ficulty. We first recalibrate this relationship in each eos- 
mological model, and then calculate the goodness of the 
relationship in that cosmology by statistics. We then 
construct a relation which is weighed by the goodness 
of each cosmology-dependent relationship, and use this 
cosmology-weighed relationship to measure the Universe. 
The procedure to calculate the probability function of a 
eosmological parameter set (denoted as f2, which includes 
both fijvf and ^Ia) is the following. 

(1) Calibrate and weigh the luminosity indicator in each 
cosmology. Given a particular set of eosmological param- 
eters (fi), we perform a multiple variable regression anal- 
ysis and get a best-fit correlation E^^iso{ri; Ep,t^). We 
evaluate the probability {w{Cl)) of using this relation as a 
cosmology-independent luminosity indicator via statis- 
tics, i.e. 



N 

E 



° 7, ISO 



(9) 



The smaller the x^, the better the fit, and hence, the 
higher probability for this cosmology-dependent relation- 
ship to serve as a eosmological independent luminosity in- 
dicator. We assume that the distribution of the xt (^) 
normal, so that the probability can be calculated as 



(10) 



(2) Regard the E^^isoi^', Ep,tfj) relationship derived in 
each cosmology as a cosmology-independent luminosity in- 
dicator without considering its systematic error, and cal- 
culate the corresponding distance modulus fi{^) [eq. 7] 
and its error ap,, which is 



x2il/2 



(12) 



(3) Calculate the theoretical distance modulus /x(0) in 
an arbitrary set of eosmological parameters (denoted by 
Q), and calculate the of /x(f2) against /t(n), i.e. 



N 



x'im-T. 



(13) 



(4) Assuming that the distribution of x^(i^|f^) is also 
normal, calculate the probability that the cosmology pa- 
rameter set O is the right one according to the luminosity 
indicator derived from the eosmological parameter set Cl, 

p{n\n) oc e->^'(^l")/2. (14) 

With eq.(10_), wc can define a cosmology- weighed likeli- 
hood by w{n)pin\n). 

(5) Integrate SI over the full cosmology parameter space 
to get the final normalized probability that the cosmology 
SI is the right one, i.e. 



p{n) = 



j^-^w{n)p{n\n)dn 



(15) 



In our calculation, the integration in eq.(15) is computed 
through summing over a wide range of the cosmology pa- 
rameter space to make the sum converge, i.e.. 



p{n) 



j:^,^w{n,)p{ni\n) 



(16) 



The essential ingredient of our method is that we do 
not include the systematical error of the Ej^iso{Sl; Ep,t^) 
relationship into cr^, ;- Instead, we evaluate the proba- 
bility that a particular relationship can be served as a 
cosmology-independent luminosity indicator using its sys- 
tematical error, and integrate over the full cosmology pa- 
rameter space to get the final probability of a cosmology 
with the parameter set fl. In Figure 4 we plot fi against //, 
with (j^, in the case of SI = 0.28 and Oa = 0.72 cosmology. 
Similar investigation could be done for other cosmologies. 
Below, we will apply the approach discussed above to in- 
vestigate the possible implications on cosmography and 
eosmological dynamics with our GRB sample. 

4.1. Implications for SIm and SI a 

In a Priedmann-Robertson- Walker (FRW) cosmology 
with mass density SIm and vacuum energy density SI a, 
the luminosity distance is given by 

Dl = c(l + z)Ho^\Slk\-^/hmn{\S}k\^^^ 

X / dz[{l + zf{l + nMz)-z{2 + z)flA]~^^^C^7) 
Jo 

where c is the speed of light, Hq is the present Hubble 
constant, Slk — ^ ~ S}m — SI a denotes the curvature of 
the universe, and "sinn" is sinh for Slk > and sin for 
Ofc < 0. For a flat universe {Slk = 0), Eq.(17) turns out to 
be c(l + z)Hq^ multiplies the integral. We calculate p{Sl,) 
with our GRB sample, where SI = {SIm,SIa)- Since both 
[(72/(1 -|- z)]^ and {(Tk/k)'^ are significantly smaller than 
the other terms in eq.(ll), they are ignored in our calcu- 
lations. Shown in Figure 5 are the most possible value of 
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(JImj^^a) and the la to Scr contours of the hkehhood in 
the (Omi^^a) plane. The most possible value of (JImj^a) 
is (0.28,0.64). The contours show that 0.05 < Qm < 0.50 
at 1(7, but Qa is poorly constrained, i.e. fl\ < 1.2 at la. 
For a flat Universe, as denoted as the dashed line in Figure 
5, the constraints are tighter, i.e. 0.13 < ftu < 0.49 and 
0.50 <Q.A< 0.85 at la. 

4.2. Implications for the cosmology dynamics 

Riess et al. (2004) found the evidence from SNe la data 
that the Universe was switched from a past decelerating 
phase to the currently accelerating phase at an epoch of 
Zt = 0.46 ± 0.13, assuming that the decelerating factor q 
evolves with redshift as q{z) = qo + zdq/dz. Following 
Riess et al. (2004), we also take q{z) = qo + zdq/dz to an- 
alyze the implications for qo and dq/dz from the current 
GRB sample. The luminosity distance in a {qo, dq/dz) 
model can be written as 



Dr = 



cjl + z) 
Ho 



(18) 



We then calculate the values of P{i^) (where = 
{qo, dq/dz)) using the cosmology- weighed standard can- 
dle method discussed above. Shown in Figure 6 arc the 
most possible values of {qo, dq/dz) and their likelihood in- 
terval contours from la to 3a. The most possible values 
of {qo, dq/dz) are (—1.0, 1.12), and at la level their values 
are constrained in the ranges of —2.23 < qo < 0.26 and 
—0.07 < dq/dz < 3.48. Although the current sample still 
does not place a tight constrain on both qo and dq/dz, it 
shows that qo tends to be less than and dq/dz tends to 
be greater than 0, suggesting that the Universe is accel- 
erating now. At a given epoch zt in the past, q{zt) = 
should be satisfied, which denotes the transition between 
the past decelerating phase and the currently accelerating 
phase. The likelihood function of Zt derived from the cur- 
rent GRB sample is shown in Figure 7. We calculate the 
best value of zt by 



Zt = 



J2pizt)zt 

Epizt) ' 



(19) 



and get = O.TSto H at la. 



5. SIMULATIONS 

We have shown that using the analysis method pro- 
posed in this paper, one can place some constraints on 
the cosmology parameters with our GRB sample. These 
constraints are, however, weaker than those obtained with 
the SNe la data, and they are of uncertainties because 
of the small GRB sample effect. To increase the signif- 
icance of the constraints, one needs a larger sample and 
smaller error bars for the measurements. In order to ac- 
cess the characteristics of the GRB sample satisfying our 
relationship observationally and how the sample size and 
the measurement precision affect the standard analysis, 
we perform some Monte Carlo simulations. We simulate a 
sample of 10"^ GRBs. Each burst is characterized by a set 
of parameters denoted as {z, Ep, S^, tb). A fluence thresh- 
old of Sth,-y = 10"'' erg s~^ is adopted. Since the observed 

^Our simulations do not sensitively depend on the Ejso distribution, 
and found that the characteristics of our simulated GRBs sample are not 



tf, is in the range of 0.4 6 days, we also require that 
ti, is in the same range to account for the selection effect 
to measure an optical lightcurve break. Our simulation 
procedures are described as follows. 

(1) Model the accumulative probability distributions of 
Ep, Eiso, and z by the observational data. We first ob- 
tain the differential distribution of these measurements. 
The Ep distribution is derived from the GRB spectral 
catalog presented by Preece et al. (2000), which is well 
modeled by dp/dlogEp cx cxp[-2(log£;p.2 - 0.38)70.452], 
where Ep^2 = Ep/lOO (Liang, Dai, & Wu 2004). The 
Eiso distribution is obtained from the current sample of 
GRBs with known redshifts. Since the Eiso distribution 
suffers observational bias at the low Eiso end, we con- 
sider only those bursts with Eiso > lO"''^"'' ergs, and get 
dp/d\ogEiso (X — 0.3 log iJjso^- The redshift distribution 
is derived by assuming that the GRB rate as a function 
of redshift is proportional to the star formation rate. The 
model SF2 from Porciani & Madau 2001 is used in this 
analysis. We truncate the redshift distribution at 10. 
Based on these differential distributions, we obtain the 
accumulative distributions, Px, where x is one of these 
parameters. We use the discrete forms of these distribu- 
tions to save the calculation time. The bin sizes of log Ep, 
logEiso, and z are taken as 0.025, 0.1, and 0.01, respec- 
tively. 

(2) Simulate a GRB. We first generate a random num- 
ber TO (0 < m < 1), and obtain the value of Xm from the 
inverse function of Px{xrn) ~ ™, i-c, Xm = ^('^)- Since 
Px is in a discrete form, we search for a bin i, which sat- 
isfies Px{xi) < TO and Px{xi+i) > to, and calculate the Xm 
value by Xm = {xi+i +Xi)/2. Repeating this step for each 
parameter, we get a simulated GRB characterized by a set 
of parameters denoted as {z, Eiso, Ep). 

(3) Calculate and examine whether or not the S-y 
satisfies our threshold setting. The gamma-ray fluence 
is calculated hy S = E[so{l + z)/AttD\{z), where Dl{z) 
is the luminosity distance at z (for a flat universe with 
r^M = 0.3). 11 S < Sth,-y, the burst is excluded. 

(4) Derive tb. We first infer a tb value from our empir- 
ical relation in a flat Universe of = 0.3, then assign 
a deviation {Atb) to the tb value. The distribution of Atb 
is taken as dN/dA\ogtb = exp{—Atl/2a), where cr = 0.1. 
This typical value is taken according to the current sample, 
which gives the mean and medium deviations as u = 0.15 
and a = 0.11, respectively. If the tb value is in the range 
of 0.4 < tb < 6 days, this burst is included in our sample. 
Otherwise, it is excluded. 

(5) Assign observational errors to Ep, S^, and tb. Since 
the observed ax/x is about 10% — 20%, we take the errors 
as ax/x = 0.25fc with a lower limit of ax/x > 5%, where 
A: is a random number between ~ 1. 

(6) Repeat steps (2) and (5) to obtain a sample of 10^ 
GRBs. 



The distributions of z, Ep, and for the simulated 
GRB sample; arc; shown in Figure 8 (solid lines). The 
observed distributions of these quantities are also shown 
for comparison (dotted lines). The observed redshift dis- 
tribution is derived from the current GRB sample with 
known redshifts (45 GRBs). The observed Ep distribution 



We have used a random distribution between 10^^'^ ~ 10^ 
significantly changed. 



ergs, 
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is taken from Preece et al. (2000). The observed S^^ is de- 
rived from the BATSE Current GRB sample'' (Cui, Liang, 
& Lu 2005). The comparisons indicate that the mock GRB 
sample tends to be a softer (low Ep) and brighter (high 
5*^) one. The redshifts of the mock GRB sample tend to 
be higher than the current GRB sample, but this might 
be due to observational biases against high redshift GRBs 
(Bloom 2003). 

We investigate the effect of the sample size on the cos- 
mological constraints with our mock GRB sample. We 
randomly select sub-samples of 25, 50, 75, and 100 GRBs 
from the mock GRB sample. We compare the la contours 
of likelihood distributions in the (Om, ^^A)-plane derived 
from these sub-samples in the left panel of Figure 9. It is 
evident that, as the sample size increases, the constraint 
on and becomes more tighter. Comparing the left 
panel of Figure 9 with the Figure 8 in Riees et al. (2004), 
we find that the likelihood contour derived from the sub- 
sample of 50 GRBs is comparable to that derived from the 
gold sample of 157 SNe la. 

Precision cosmology requires accurate observations. 
Modern sophisticated observation techniques in distant 
SNe la (e.g. Riess et al. 1998; Schmidt et al. 1998; Perl- 
mutter et al. 1999) and cosmic microwave background 
(CMB) fluctuations (e.g. Bennett et al. 2003; Spergel 
et al. 2003) have made great progress on modern preci- 
sion cosmology. We inspect the uncertainties of the dis- 
tance modulus derived from the SNe la data, and find 
that the average uncertainty is a-jjM ~ 0.25, while for our 
GRB sample it is 0.45. Increasing observational precision 
(i.e. reducing the errors) should significantly improve the 
constraints on the cosmological parameters. We simulate 
another GRB sample with systematically smaller obser- 
vational errors, i.e. a^/x = 0.15A; in the step 5 of our 
simulation procedure. We get a sample with aoM ^ 0.28, 
comparable to the SNe la gold sample. The comparison of 
the likelihood contours (la) in the D,m — f^A plane derived 
from a sample of 50 mock GRBs with aoM ^ 0.45 (the 
line contour) and with aoM ~ 0.28 (the grey contour) are 
shown in the right panel of Figure 9. It is found that the 
latter is significantly tighter, which is comparable to that 
derived from a sample of 100 mock GRBs with an average 
error in modulus of 0.45. 

The results in Figures 9 indicate that tighter constraints 
on cosmological parameters could be achieved by either 
enlarging the sample size or increasing the observational 
precision. If a sample of 50 GRBs with a comparable ob- 
servational precision as SNe la gold sample could be estab- 
lished, the constraints are even tighter than those derived 
from the SNe la gold sample. 

6. CONCLUSIONS AND DISCUSSION 

Without imposing any theoretical models and assump- 
tions, we investigate the relationship among E^^iso, E^ 
and using a multiple variable regression method. Our 

GRB sample includes 15 bursts, whose E^ and are well 
measured. The results indicate that £^7,130 strongly de- 
pends on both Ep and t^^ with a very small dispersion, e.g. 
eq.(5) for a flat Universe with Hm = 0.28. We also per- 
form a similar analysis by replacing E^^iso by the isotropic 

^http://cossc.gsfc. nasa. gov /batse/ 



afterglow energies in the X-ray and optical bands, and 
find that these energies are essentially not related to E^ 

and at all. We then use the E^,isoiEp,t-^^) relationship 
as a luminosity indicator to infer the possible cosmology 
implications from the GRB sample. Since this relation- 
ship is cosmology-dependent, we suggest a new method to 
weigh various cosmology-dependent relationships with its 
probability of being the right one, and use the cosmology- 
weighed standard candle to explore the most plausible cos- 
mological parameters. Our results show that the most pos- 
sible values are (Qmj^a) = (0.28,0.64). At la level, we 
have 0.05 < J^m < 0.50 and J^a < 1-2. In the case of a 
flat Universe, the la constraints are 0.13 < flu < 0.49 
and 0.50 < f^A < 0.85. The decelerating factor of the Uni- 
verse (q) and its cosmological evolution (dq/dz) are also 
investigated with an evolutionary form oi q = qo + zdq/dz. 
The GRB sample implies that the most possible values of 
{qo, dq/dz) are (—1.00,1.12), and they are constrained in 
the ranges of -2.23 < go < 0.26 and -0.07 < dq/dz < 3.48 
at la level. A transition redshift between the deceleration 
and the acceleration phases of the Universe is inferred as 
Zt = 0.78to;23 at la level from the GRB sample. 

As a Imninosity indicator, our model-independent 
Ej^iso{Ep,t^) relationship takes the advantage upon the 
previous Ghirlanda-relation in that only pure observa- 
tional data are involved. Since this luminosity indicator 
is cosmology-dependent, we use a strategy through weigh- 
ing this relationship in all possible cosmologies to statis- 
tically study the cosmography and cosmological dynam- 
ics. A similar method has been used in the SN cosmology 
when dealing with the uncertainty of the present Hubble 
constant Hq. In their method (e.g. Riess et al. 1998), the 
systematic error of Hq is not included to calculate the error 
of the distance modiflus. Rather, they integrated the prob- 
ability of Hq over a large range values (without weighing 
for each value of Hq). This is the so-called marginalization 
method. We also perform this marginalization method to 
deal with our coefficients (kq, ki and K2), and re-do the 
cosmology-analysis. This is equivalent to integrating over 
the whole cosmology parameter space without weighing, 
i.e., 

p{n) = / p{n\n)dQ. (20) 

Jn 

The result using this method to constrain Om and CIa is 
presented in Figure 10. Comparing it with Fig. 5, we find 

that both methods give consistent results, but Fig. 5 gives 
a tighter constraint on cosmological parameters. This is 
understandable, since the weighing method reduces the 
contributions of side lobe around the "true" cosmologies. 
In any case, an essential ingredient of both methods is that 
the uncertainty of the standard candle itself is not included 
in calculating the uncertainty of the distance modulus de- 
rived from the data. If the uncertainty of the standard 
candle is indeed included in the uncertainty of the distance 
modulus, with eqs.(15) and (20) to calculate p(0), one gets 
a very loose constraint (Fig. 11). Even at la level the cur- 
rent GRB sample cannot place any meaningful constraints 
on both flu and ^Ia)- We believe, however, that in such a 
treatment, the uncertainty of the distance modiflus is over- 
estimated, since the error introduced from measurements 
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should not be mixed with the systematic uncertainty of 
the standard candle. 

The GRB sample from which our relationship is drawn 
is currently small. The constraints on the cosmological pa- 
rameters derived from this sample arc weaker than those 
from the SNe la gold sample. Our simulations indicate 
that either the enlargement of the sample size or the in- 
crease of the observational precision could greatly improve 
the constraints on the cosmological parameters. A sample 
of 50 bursts with the current observational precision would 
be comparable to the 157 SNe la gold sample in constrain- 
ing cosmology, and a better constraint is achievable with 
better observational precisions or an even larger sample 
size. 

Our simulations also indicate that the GRB sample sat- 
isfying our relationship observationally tends to be a soft 
and bright one, for which tj, is in the reasonable range for 
detection. Detailed optical afterglow light curves cover- 
ing from a few hours to about ten days after the burst 
trigger"* are required to measure the tt value. The ob- 
served tb ranges from 0.4 days to 5 days in the current GRB 
sample. In the CGRO/BATSE duration table^, there are 
- 1500 long GRBs. To test the probabihty of a BATSE 
burst having a th in the range of 0.4 days to 5 days, wo 
perform a simulation similar to that described in §5, but 
take the Ep and Sj distributions directly from the BATSE 
observations. We find that the probability is ~ 0.15 in the 
cosmology of Hm — 0.3 and n\ — 0.7. Among well local- 
ized GRBs about 50% of bursts are optically bright. We 
thus estimate that there have been ~ 110 BATSE GRBs 
that might have been detected to satisfy our E~,_[so(Ep,ti^) 
relationship. As shown in Figure 9, such a sample is com- 
parable to the SNe la gold sample for constraining cos- 
mological parameters. A dedicated GRB mission carrying 
a BATSE-liked GRB detector and having the capability 
of precisely localizing and following up GRBs (like Swift) 
would be ideal to establish a homogenous GRB sample to 
perform precision GRB cosmology (Lamb et al. 2005b). 

Since launched on Nov. 20, 2004, the Swift mission 
(Gehrels et al. 2004) is regularly detecting GRBs with 
a rate of ^ 80 bursts per year. Detailed X-ray and 
UV/optical afterglow observations spanning from 1 minute 
to several days after the burst have been performed for 
most of the bursts. However, the energy band of the Swift 
Burst Alert Telescope (BAT) is narrow, i.e. 15-150 keV. 
As we shown in Figure 8, the typical Ep of a burst in 
our sample is marginally at the end of the BAT energy 
band. As a result, BAT is not ideal for the purpose of 
expanding the sample for GRB cosmology. As a result, 
we do not expect a dramatic enlargement of our sample 
in the Swift era. Nonetheless, those bursts with an Ep of 
~ (50 — 100) keV^ could have their Ep well-measured by 
Swift. We therefore highly recommend a detailed optical 
follow-up observations for these bursts with UVOT and/or 
other ground-based optical telescopes. This would present 



an opportunity to enlarge our sample with the Swift data. 

The major advantage of GRBs serving as a standard 
candle over SNe la is their high-redshift nature. The ob- 
served spectra and fluences of high redshift GRBs may 
not be too different from the nearby ones (Lamb & Re- 
ichart 2000). For example, the fluence of GRB 000131 
(z = 4.5) is 1 X lO"-"^ erg cni"^ ^^le 25-100 keV band 
(Hurley et al. 2000), which is significantly larger than 
the fluence of typical GRBs (~ 10~^ erg cm~^ in the 25- 
2000 keV band). The highest redshift burst in our sam- 
ple is GRB 020124 (z 3.2). Its fluence is 6.8 x 10"^ 
erg cm~^ in the 30-400 keV band, and its observed peak 
energy is 120 keV. These indicate that the current GRB 
missions, such as Swift and HETE-2, are adequate to 
observe high redshift GRBs^. We explore how the con- 
straints on cosmological parameters are improved by iden- 
tifying several high-redshift bursts. We artificially select 

5 high redshift GRBs from our simulated GRB sample 
with z ^ 4.0, 5.0, 6.0, 7.0 8.0, respectively, and with ob- 
servational errors a^/x ~ 0.25fc. The constraints on the 
cosmological parameters from these pseudo high-z GRBs 
together with the current observed GRB sample are shown 
in Figure 12 (the grey contour), where the results from the 
current observed GRB sample are also plotted for com- 
parison (the line contours). It is found that adding a few 
high-z GRBs could result in much tighter constraints on 
cosmological parameters. Identifying high-z GRBs and 
measuring their Ep and f'j, are therefore essential for the 
GRB cosmology in the near future. 

Our model-independent relationship is close to the 
Ghirlanda relationship, which was derived based on a sim- 
plest version of GRB jet models. In such a model invoking 
a jet with energy uniformly distributed in the jet cone, 
the observable tb is physically related to the epoch when 
the bulk Lorentz factor of the ejecta is reduced to the in- 
verse of the jet opening angle (Rhoads 1999; Sari et al. 
1999). Relating tb to the jet openning angle, the jet en- 
ergy is then given by £jct oc {Ej^isot^)'^^^{nriy y/^, where n 
is circum-burst medium density and r]^ is the efficiency of 
GRBs. The Ghirlanda-relation can be then expressed as 
ii^iso fx Ept^^{nr]^)~^^^ . Comparing this with our model- 
independent relationship (eq.5), we can see that both rela- 
tions are roughly consistent with each other if n and r]j are 
universal among bursts. As discussed above, the motiva- 
tions for us to introduce our multi- variable relationship are 
two folds. Firstly, n and r]^ may not be constant and actu- 
ally vary from burst to burst. This introduces a lot more 
uncertainties in the Ghirlanda relationship (e.g. Friedman 

6 Bloom 2005). Secondly and more importantly, there 
is no straightforward interpretation of the relation within 
the jet model. Jumping out from the jet model framework 
would give more freedom of theoretical interpretations. 

The tight relation of Eiso{Ep, t^) is very intriguing, and 
its physical reason calls for investigation. The fact that 



"^Starting from about ten days, the contributions from the underlying SN and host galaxy components may become prominent, and the 
afterglow level may be too faint to be detected. 

® http://gammaray.msfc.nasa.gov/batse/grb/catalog/current 
®Such a burst tend to be an X-ray rich GRB (Lamb et al. 2005a). 

^Strictly speaking, wc refer to the optical band in the cosmological rest frame to define This is not an issue if the GRB redshift is 
small. For high-2; GRBs, the optical band in the observer's frame is highly extincted by neutral hydrogen, but one could still detect tb from 
the infrared band. Infrared band observations are also essential for identify high-z GRBs. IR follow up observations are therefore essential to 
add the high-z bursts in our sample. 
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Ej^iso strongly depends on and t^, while both -ExA,iso 

and -EoAjiso do not, implies that is a quantity related to 
GRBs rather than to their afterglows. A similar signature 
was previously found by Salmonson & Galama (2002), who 
discovered a tight correlation between the pulse spectral 
lag of GRB light curves and tb. We therefore suspect that 
tb might be a unique probe for the GRB prompt emission 
properties. Within the jet scenario, the anti-correlation 
between tb and i^^^iso (first revealed by Frail et al. 2001) 
may be physically related to the different metalicity abun- 
dances of the progenitor stars (e.g. metal-poor stars ro- 
tate more rapidly, and the GRBs they produce are more 
energetic and have more coUimatcd jets, MacFadycn & 
Woosley 1999; Ramirez-Ruiz et al. 2002) or may be simply 
a manifestation of the viewing angle effect in a structured- 
jet scenario (Rossi et al. 2002; Zhang & Meszaros 2002b). 
Such an anti-correlation, when combined with the physical 
models of i?-y.iso — Ep correlations (e.g. Zhang & Meszaros 
2002a; Rees & Meszaros 2005), may be able to inter- 



pret the observed Eiso{Ep, t^,) relation, although a detailed 
model is yet constructed. Alternatively, there might be a 
completely different physical reason under the Eiso{Ep, t^) 
relationship which is not attached to the jet picture. One 
possibility is that the spectral break in the prompt emis- 
sion and the temporal break in the optical band may be 
related to a same evolving break in the electron spectral 
distribution (B. Zhang, 2005, in preparation). In such 
an interpretation, the temporal break time in the optical 
band is expected to be different from those in the radio or 
in the X-ray bands. Since so far there is no solid proof for 
the achromatic nature in broad bands for any "jet break" , 
such a possibility is not ruled out. 

We are grateful to the anonymous referee for valuable 
comments. Wc also thank Z. G. Dai, D. Xu, Y. Z. Fan, 
X. F. Wu, G. Rhee, D. Lamb, G. Ghirlanda, D. Lazzati, 
T. Piran, Y. P. Qin, and B. B. Zhang for helpful discus- 
sions. This work is supported by NASA NNG04GD51G, 
a NASA Swift GI (Cycle 1) program, and the National 
Natural Science Foundation of China (No. 10463001). 
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Table 1 

Prompt emission parameters of the GRB sample adopted in this paper 



GRB 


z 


(3)(keV) 


a 


/? 


P=Y / \ 

Ot(<ts) 


Band 


Rcfs. 


(1) 


(2) 


(4) 


(5) 


(6)(erg.cm~^) 


(7)(keV) 


(8) 


980703 


0.966 


254(50.8) 


-1.31 


-2.40 


22.6(2.3) 


20-2000 


1; 


2; 2; 2 


990123 


1.6 


780.8(61.9) 


-0.89 


-2.45 


300(40) 


40-700 


3; 


4 ; 4;4 


990510 


1.62 


161.5(16.1) 


-1.23 


-2.70 


19(2) 


40-700 


5: 


4; 4; 4 


990712 


0.43 


65(11) 


-1.88 


-2.48 


6.5(0.3) 


40-700 


5: 


4: 4: 4 


991216 


1.02 


317.3(63.4) 


-1.23 


-2.18 


194(19) 


20-2000 


6; 2; 2; 2 


011211 


2.14 


59.2(7.6) 


-0.84 


-2.30 


5.0(0.5) 


40-700 


7; 8; 8; 7 


020124 


3.2 


86.9(15.0) 


-0.79 


-2.30 


8.1(0.8) 


2-400 


9: 10: 10: 10 


020405 


0.69 


192.5(53.8) 


0.00 


-1.87 


74.0(0.7) 


15-2000 


11 


11: 11: 11 


020813 


1.25 


142(13) 


-0.94 


-1.57 


97.9(10) 


2-400 


12 


10: 10; 10 


021004 


2.332 


79.8(30) 


-1.01 


-2.30 


2.6(0.6) 


2-400 


13 


10; 10; 10 


021211 


1.006 


46.8(5.5) 


-0.86 


-2.18 


3.5(0.1) 


2-400 


14 


10; 10; 10 


030226 


1.986 


97(20) 


-0.89 


-2.30 


5.61(0.65) 


2-400 


15 


10: 10: 10 


030328 


1.52 


126.3(13.5) 


-1.14 


-2.09 


37.0(1.4) 


2-400 


16 


10; 10; 10 


030329 


0.1685 


67.9(2.2) 


-1.26 


-2.28 


163(10) 


2-400 


17 


10; 10; 10 


030429 


2.6564 


35(9) 


-1.12 


-2.30 


0.85(0.14) 


2-400 


18 


10; 10; 10 



References. — References are in order for z, E'^^ , [a,/?], S'7:(l) Djorgovski et al. 1998; (2) Jimenez et aL 2001; (3) Kulkarni et al. 1999; 
(4) Ainati et al. 2002; (5) Vreeswijk et al. 2001; (6)Djorgovski et al. 1999; (7) Holland et al. 2002; (8) Amati 2004; (9) Hjorth et al. 2003; 
(lO)Sakainoto et al. 2004b; (11) Price et al. 2003; (12) Bartli et al. 2003; (13) MoUer et al. 2002; (14) Vreeswijk et al. 2003; (15) Greiner et 
al. 2003; (16) Martini et al. 2003; (17)Bloom et al. 2003c; (18) Weidinger et al. 2003. 



Table 2 

x-ray and optical afterglow data of the grb sample adopted in this paper 



GRB 




Epoch^ 






ib(atj'' 
(6) (days) 






Ref.*^ 




(1) 


(2) 


(3) (hours) 


(4) 


(5) 


(7) 


(8) 


(9) 


(10) 


980703 


1.24 


34 


4 


18.24(4.97) 


3.4(0.5) 






1 


20.1 


990123 


1.08 


6 


110 


66.09(6.33) 


2.04(0.46) 


1.17 


1.57 


2 


19.4 


990510 


1.41 


8.7 


47.8 


41.07(3.68) 


1.6(0.2) 


0.46 


1.85 


3 


18.1 


990712 










1.6(0.2) 


0.83 


3.06 


4 


19.5 


991216 


1.61 


4 


1240 


287.21(14.73) 


1.2(0.4) 


1 


1.8 


5 


16.9 


011211 


1.5 


11 


1.9 


2.23(0.39) 


1.56(0.02) 


0.95 


2.11 


6 


20.1 


020124 










3(0.4) 






7 


21.6 


020405 


1.15 


41 


13.6 


68.98(20.21) 


1.67(0.52) 


1.4 


1.95 


8 


18.3 


020813 


1.42 


39 


22 


113.98(17.01) 


0.43(0.06) 


0.76 


1.46 


9 


19.1 


021004 


1.56 


20.81 


4.3 


13.5(2.47) 


4.74(0.14) 


0.85 


1.43 


10 


18.4 


021211 










1.4(0.5) 






11 


21.3 


030226* 




37.1 


0.32 


12.3 


1.04(0.12) 


0.77 


1.99 


12 


19.5 


030328* 




15.33 


3 




0.8(0.1) 


1.0 


1.6 


13 


20.2 


030329® 


1.74 


4.85 


1400 


467(23) 


0.5(0.1) 


1.18 


1.81 


14 


14.7 


030429 










1.77(1) 


0.88 


2.87 


15 


20.2 



"Temporal decay index and X-ray afterglow flux in 2-10 keV band at a given observed epoch, i^j^.ioh is the extrapolated/interpolated X-ray 
afterglow flux at 10 hours after the GRB trigger. The fluxes are in units of lO^^"* ergs cni^^ s They are taken from Berger et al. (2003) 
except for those with marks: 030226 (Pedersenct al. 2003); 030328 (Butler et al. 2003); 030329 (Marshall & Swank 2003; Marshall, Markwardt, 
& Swank 2003; Tiengo, Mereghetti, & SchartelA 2003a, b) 

''Temporal break (error) and temporal indices before and after the break, and their references: (1) Frail et al. 2003; (2) Kulkarni et al. 1999; 
(3) Stanek et al. 1999; (4) Bjornsson et al. 2001; (5) Halpern et al. 2000; (6) .lakobsson et al. 2003; (7) Berger et al. 2002; (8) Price et al. 
2003; (9) Barth et al. 2003; (10) Holland et al. 2003; (11) Holland et al. 2004; (12) Klose et al. 2004; (13) Andersen et al. 2003; (14) Berger 
et al. 2003; (15) Jakobsson et al. 2004a. 

'^R-band magnitude adjusted to 11 hours after the burst trigger (from Jakobsson et al. 2004b). 
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Table 3 

Derived isotropic energies, rest frame peak energies and rest frame temporal breaks for the GRB sample 

ADOPTED IN THIS PAPER (ASSUMING Qm = 0.28 AND f^A = 0.72) 



GRB 


(l)(crg) 


logi;xA,iso 

(2) (erg) 


logiJoA.iso 

(3)(erg) 


(4)(kcV) 


1 ,' / \ 
(5) (day) 


980703 


52.85(0.04) 


47.69 




2.70(0.09) 


0.238(0.064) 


990123 


54.64(0.06) 


48.80 


46.09 


3.31(0.03) 


-0.105(0.098) 


990510 


53.29(0.05) 


48.58 


46.60 


2.63(0.04) 


-0.222(0.008) 


990712 


51.88(0.02) 




44.08 


1.97(0.07) 


0.049(0.054) 


991216 


53.85(0.04) 


49.00 


45.70 


2.81(0.09) 


-0.226(0.145) 


011211 


53.01(0.04) 


47.72 


45.84 


2.27(0.06) 


-0.304(0.006) 


020124 


53.37(0.05) 






2.70(0.08) 


-0.146(0.058) 


020405 


53.17(0.01) 


47.92 


45.91 


2.51(0.12) 


-0.005(0.135) 


020813 


54.13(0.06) 


48.90 


45.42 


2.68(0.09) 


-0.719(0.061) 


021004 


52.66(0.10) 


48.54 


46.82 


2.42(0.16) 


0.153(0.013) 


021211 


52.05(0.03) 






1.97(0.05) 


-0.156(0.155) 


030226 


52.90(0.05) 


47.52 


46.14 


2.46(0.09) 


-0.458(0.050) 


030328 


53.60(0.02) 


47.66 


45.51 


2.50(0.05) 


-0.498(0.054) 


030329 


52.19(0.04) 


48.27 


45.40 


1.90(0.01) 


-0.369(0.087) 


030429 


52.24(0.07) 




46.30 


2.11(0.11) 


-0.315(0.245) 



Table 4 

The RESULTS of multiple variable regression ANALYSIS (Assuming Om = 0.28 and = 0.72) 







-E'x,iso(-Ep,tf,) 


£'R,iso(-Ep,tj,) 




48.0±0.4(< 10^^) 


46.27 ± 1.35(< lO-'') 


44.22 ± 1.42(< 10"*) 




1.94±0.17(< 10-"^) 


0.74 ±0.51(0.18) 


0.64 ±0.57(0.26) 


K2(Pt) 


-1.24 ±0.23(2 X 10"^) 


-0.30 ±0.62(0.64) 


0.29 ±0.88(0.75) 


Global F-tcst statistics 


115.4 


1.08 


0.77 


and probability pp 


< 10-* 


0.39 


0.49 


Global correlation r 


0.96 ±0.21 


0.46 ± 0.24 


0.38 ± 0.28 


and probability Ps 


< 10-4 


0.15 


0.22 
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Fig. 1. — Log E-f^iao calculated by the empirical relationship from our multiple variable regression analysis as compajred with log E-^^iso 
derived from the observed fluence with the cosmological parameters of Qm = 0.28 and Qa = 0.72. The solid line is the regression line for the 
two quantities. 
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Fig. 2. - Log i/'xA.iso calculated by the empirical relationship from our multiple variable regression analysis as compared with log BxA,ii 
derived with the cosmological parameters of VIm = 0.28 and Oa = 0.72. 
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Fig. 3. — Log i?oA,iso calculated by the empirical relationship from our multiple variable regression analysis as compared with log i7oA,iso 
derived with the cosmological parameters of CIm = 0.28 and Ha = 0.72. 
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Flo. 4. — The distance modulus derived from tiie data, fi, and its observational error, ct^, are plotted against the distance modulus derived 
from theory, jj,. The cosmological parameters are adopted as Qm = 0.28 and Qa = 0.72. 
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Fig. 5. — The contours of likelihood interval distributions in the (Qm, f2A)-plane inferred from the current GRB sample using the method 
developed in §4. The cross marks the most possible value of {Qm,^a), which is (0.28,0.64). The contours give 0.05 < Qm < 0.50(1it). 
Considering a flat universe (the dashed line), the contours yield 0.13 < flM < 0.49 and 0.50 < fl\ < 0.85 (1<t). 
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Fig. 6. — The contours of likelihood interval distributions in the (qo, c(<;/(i2)-plane inferred from the current GRB sample using the method 
developed in §4. The most possible values of {qQ,dq/dz) are (—1.00,1.12) (the cross). At la level their values are constrained in the ranges 
of -2.23 < go < 0.26 and -0.07 < dq/dz < 3.48, respectively. 
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Fig. 8. — The distributions of z, Ep, and Sj for the simulated GRB sample satisfying our model-independent standard candle relationship. 
For comparison, the imposed dotted lines are the distributions derived from the observational data. 



19 




Fig. 9. — A comparison of the la likelihood contours for different simulated samples: right — Simulations for different sample sizes: dotted 
contour — 25 GRBs, solid contour — 50 GRBs, dashed contour — 75 GRBs, grey contour — 100 GRBs. The same observational errors 
{<Tx/x = 0.25fc, k is a random number between ~ 1) are adopted; left: Simulations for a same sample size (50 GRBs) but for different 
observational errors: the line contour — Tx/x = 0.25fc and the grey contour — rix/x = 0.15fe (see the procedure of our simulations). The 
dotted line is for a flat Universe. 
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Fig. 10. — The contours of likelihood interval distributions in the {CIm, nA)-plane derived by the marginalization method. The cross marks 
the most possible value of (f2M, S^a) , which is (0.28, 0.52). The contours give 0.05 < Qm < 0.61(1it). Considering a flat universe (the dashed 
Une), the contours yield 0.14 < Hm < 0.58 and 0.40 < f^A < 0.84 (la). 
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Fig. 11. — Same as Fig. 5 but the uncertainties of tlie parameters in the empirical relationships are also included into the error of the 
distance modulus in the calculation of p{il) by Eq. 15. Only Icr interval is shown. 
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Fig. 12. — A demonstration of the potential constraints on cosmological parameters with high redshift GRBs. The gray contours are the 
results derived from 5 pseudo high redshift GRBs together with the observed ORB sample [Itr (dark gray) to 3a (light gray)], and the line 
contours are the results from the current observed GRB sample only. 



